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Abstract. The interaction of partially ionized plasmas with an electromagnetic field is 
investigated using quantum statistical methods. A general statistical expression for the current 
density of a plasma in an electromagnetic field is presented and considered in the high field 
regime. Expressions for the coUisional absorption are derived and discussed. Further, partially 
ionized plasmas are considered. Plasma Bloch equations for the description of bound-free 
transitions are given and the absorption coefficient as well as rate coefficients for multiphoton 
ionization are derived and numerical results are presented. 



1. Introduction 

Due to the enormous progress in short-pulse laser physics, the interaction of matter with 
electromagnetic fields is a topic of rapidly growing interest [1, 2]. In this paper we consider 
a partially ionized plasma under the influence of an external laser field. 

A partially ionized plasma is characterized by the densities of its constituents: electrons 
(ne), ions (n^), and atoms (ra^). The temporal change of these densities due to ionization and 
recombination is determined by rate equations. The rate equation for the atom density has the 
following form 



Here, a and /3 are the collisional rate coefficients for the reaction A + e ^ e + e + i, and a and 
are the radiative rate coefficients according to A + fko ^ e + i. In the field-free case, and 
P^ vanish. The collisional coefficients for dense plasmas, intensively discussed in [3, 4, 5], are 
essentially determined by the interaction. The density dependence of a can be seen in Fig. 1 



In the stationary case, from the rate equation a mass action law follows which determines 
the equilibrium composition of the partially ionized plasma. This is shown in Fig. 2 [3]. A 
characteristic feature is the strong increase of the degree of ionization (oe), well-known as density 
ionization (Mott effect). 

Under the action of an external electromagnetic field given by, e.g.. 



duA 



PneUeni — aUeUA + P UeTli — a UA- 




dt 



[3]. 



A{t) = Aosin(a;t), E(t) = -— A(t) = Eocos(u;t), Eq = -wAq, 



(2) 




Figure 1. Impact ionization coefficient for the atomic 
ground state in a hydrogen plasma as a function of the 
free electron density at several temperatures. Dashed 
lines: simple analytical formula for the ionization cross 
section proposed in [4], solid lines: using numerical results 
for the ionization cross section [5, 3]. 



the situation is more complicated. Several additional interesting processes occur, e.g. (i) 
field induced transitions between bound and scattering states, i.e. radiation ionization and 
recombination described by and (3^, (ii) excitation and deexcitation processes of atoms in 
the plasma, and (iii) processes between scattering states like (inverse) bremsstrahlung (collisional 
absorption) . 

All these processes are connected with an energy exchange between plasma and field which 
is determined by the electrical current density j, 



dt 



dt 



(3) 



i.e., the change of the total energy of the system of particles is equal to j ■ E which is in turn the 

energy loss of the electromagnetic field due to Poynting's theorem. From the quantity j wc get, 
as is well known, further important quantities like the polarization dP /dt = j, the absorption, 
reflection and refraction coefficients. 
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Figure 2. Equilibrium composition of a hydrogen plasma at T = 15000 K where aa denotes 
the fraction of the density of free particles of the species a to the total electron density. 



A central quantity, therefore, is the electrical current density defined by 



(4) 



where fa{Pat) is the single-particle gauge invariant Wigner function of the species a. For the 
determination of /„ we start from the kinetic equation [6, 7] 

d 1 

- + e„E(i) • Vfc„ } fa{K, t) = laiK, t) , (5) 



at 

where the collision integral is given by (V being the volume and Vab the interaction potential) 

Ia{Pa,t) = X^yy (^27rh)?> iPaPb\[Vab, Fab{t)]\PbPa) 

S-^ I f (iq(ipi(ip2, , . , 

= L V y (2vrn)6 t'^^P" - Pi - q) - '^(P« - Pi)] 

xK6(q)i^o6(pi,Pi + q,P2,P2 - q,*) • (6) 

The two-particle density matrix is connected with a two-particle correlation function by 

Fabit) = {ihf9<,im, 1'2't) = {ihfg< {It, I't) g< {2t, 2't) ± exch. + {ihfg<''^^^{12t, l'2't) , 

where we denoted 1 = ri, ,si etc. For special physical situations, it is easier to derive appropriate 
approximations for other functions. In place of the correlation part of g^^, we can put any single- 
time two-particle correlation function, especially also L^j^ — which is the correlation function 
of density fluctuations [8, 9], 

(i^)L<(ll't,22Y) = {Sphi22't')5pa{n't)) 

with Spai^l't) = ^+{l',t)^ail,t)-{^t{^',t)^aM) (7) 

In the momentum balance, there are contributions from collisions between different species only. 
For the case a ^ b, we have (a = e, i) 

Fabit) = iinfg<ilt,rt)g<{2t,2't) + iih)L<{n't,2'2't). (8) 

Explicit expressions for the collision integral were given in the classical case first by Silin 
[10] and later by Klimontovich [11] and in the quantum mechanical case in the papers [6, 12]. 
In generalization to the well-known Boltzmann-like collision integrals, they describe short-time 
processes, i.e., they are non-Markovian [13, 14, 15, 16]. 

In a partially ionized plasma, we have to account for bound and scattering states and 
transitions between them caused by the external field. Therefore, it is useful to expand Fat with 
respect to the eigcnstatcs of the unperturbed two-particle Hamiltonian (the spin dependence is 
not written explicitly) 

i^a6(pi,Pi + q,P2,P2-q,0 = Yl (PiP2|«P)<Pa|i^a6(i)|a'P'>(P'«'|Pi + q,P2-q) (9) 

Pa,F'a' 

with a = {n,l,m} for bound states and a = p for scattering states. P is the center-of- 
mass momentum. As we will show lateron, due to the external field, there are nonvanishing 
nondiagonal elements F^" . They are connected with transitions between different states a and 
a'. 

The balance equation for the electrical current density follows in well-known manner from 
Eq. (5) [8] 



Further progress is connected with the determination of the functions or F^b, respectively. 

In the following, we will consider the free-free transition using approximations for in the 
next section. Then the current is calculated from Eq. (10), and the collisional absorption (inverse 
bremsstrahlung) in the linear response regime as well as in the high-field case is determined and 
discussed. 

In the third section, bound states arc included. The equation of motion for is the Bethe- 
Salpeter equation (BSE) for plasmas which has been discussed in detail in [17]. We will derive 
plasma Bloch equations which form a system of equations for the determination of the matrix 
elements of (9). Here, the equations for the diagonal elements describe the time evolution of the 
occupation of bound and scattering states, respectively, driven by the non-diagonal elements. 
The coupled system, therefore, describes the ionization kinetics of a partially ionized plasma. 



2. Collisional absorption in a fully ionized plasma 

Let us first consider the energy exchange between field and plasma due to transitions between 
scattering states. This problem has been dealt with in several papers, first for classical plasmas, 
e.g. [18, 10, 11, 19, 20]. Quantum mechanical treatments were given in more recent papers, for 
the strong field case, e.g., [21, 22, 23], and in linear response theory [24]. 

Starting point for our investigation is Eq. (10). The r.h.s. of this equation is determined 
by the function L^. This quantity follows from the Bethe-Salpeter equation on the Keldysh 
contour 

L„(.(12,l'2') = n„,,(12,l'2') + ^ / d3d4nae(13,l'3+)Kd(34)Ld;,(42,4+2'). (11) 
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In a plasma in a strong laser field, an approximation in lowest order of Vie is appropriate because 
the coupling between species with different charges is weaker than their coupling to the field. 
The coupling between particles with equal charges in the subsystem, however, is not affected 
by the field. Then the so-called polarization functions Ilab can adopted to be approximately 
diagonal, Ilab = ^ab^a- The solution has now the following structure (for brevity all arguments 
are suppressed) 
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(12) 



Going back from the contour to the physical time axis using the Langreth-Wilkins rules we find 
for the correlation functions of the density fluctuation [9] 

Lj.(q; t, t') = j di [£|(q; t, f) Vei{q) Al(q; t, t') + £f,(q; t, t) Vei{q) 4(q; *')] ■ (13) 



For the electron current (10) there follows that 
.2 r ^3 



%Ut)-ne^^{t) = Re 
dt me 



l\z -^K»(g)27ri / dt 

{2-KhY meh Jto 



See{ci;t,t)Vei{q)ji:ti{ci;t,t) 



+C^e{^;t,t)Vei{q)Su{ci;t,t) 

where we introduced the dynamical structure factor 

ih r 

27r5a„(q; t,t) = - [/:> (q; t, t) + (q; t, t) 



(14) 



(15) 



The functions in the collision term depend on the current and the electrical field, respectively. 
This dependence can be made explicit if one assumes that each subsystem (electrons and ions) 



is in local equilibrium with a temperature with respect to a coordinate frame moving with 

the mean velocity Ua{t) [3, 25]. The transformation between such a coordinate system and a 
system at rest is given by f = r — J^* dtUa{t). The Fourier transforms in the two coordinate 
systems are connected by 

>C„„(q,M2) = e-i'i-^*2^*"<^(*")4„(q,ti-t2), (16) 

where Caa denotes the local equilibrium function depending on the time difference only. One 
gets (omitting the tilde from now on) 

d f dPa e 27r /"* _ r 

-i{t) = e^u:lm-'^^]j^^:;^yie{Q)-j-j^di[^^^^ (17) 

x£|t(q;f-t) + £f,(q;t-f)yei(9)<5u(q;i-i)lexp|-^^q- / dtu{h)\ ■ 

J fin^e^ j-f. J 

Thus the source term in the current balance equation depends on the current itself in a non-linear 
way. 

We will consider this non-linear equation in the following in two limiting cases. 
2.1. Linear response case 

Let us first consider the energy exchange in a weak laser field. To describe the situation we 
introduce the quiver velocity 

vo = . (18) 

Then the weak field case is characterized by vq /uth <C 1 where vth is the electron thermal velocity. 
Therefore, the situation corresponds to the linear response regime. In this case we can expand 
the exponential function in Eq. (17), 

expj-^^q- / dtij(ii)| Ril--l^q. / dtii{h)- (19) 

Adopting a harmonic time dependence of the electric field, we have for the current j{t) = 
j(a;)e~"^* -|- j*(a;)e*'^* and get from the above equation 

jM = . """^^ , , E(a;), (20) 
which is a generalized Drude equation with a complex electron-ion collision frequency Ugi given 

by 

oo 

uU^)=i - / dqq'Su{q)V{q)\e;,\<i;u^)-s;,\<i;0)] (21) 



with the dielectric function of the electron subsystem 

e-i(q;a;) = I + e^V{q) Cee{ci,u:) , (22) 

the static ion- ion structure factor Sii{q), . Using for the RPA, this is the well-known 
expression first given by Bekefi [26] and later for instance by Reinholz et al. [24]. One gets 
immediately the dynamical conductivity 



Figure 3. Electron-ion collision frequency as a 

function of the laser frequency for a hydrogen 
plasma (ng = 10^^ cm'^; T = 10^ K) for 
two different field strengths. The upper curve 
of each pair corresponds to the full dynamical 
screening, the lower one to the static screening 
approximation. 

lo" lo' 

In the limit of high frequencies, one gets cfiuj) k, [eQio^/io^)[iu + Vei{u)\ [18], and the averaged 
absorbed energy is given by 

(j • E) = 1 / dt'iit') ■ E(i') = ^Re ^ = ^Re (eoE^) . (24) 

The dependence of the quantity Re i^gj on the laser frequency is shown in Fig. 3 for two different 
relations ^'o/^^th• Only the upper set of curves correspond to the linear response formula (21). 
Obviously, statical and dynamical screening exhibit a qualitatively different behavior. The 
dynamical screening leads to a maximum in the vicinity of the plasma frequency, u ujp. The 
lower set of curves is beyond the linear response regime. This strong field case will be considered 
in the next paragraph. 

2.2. Strong high-frequency fields 

The strong field case is characterized by fo/^^th ^ 1- For the first time such a situation was 
discussed by Silin [10] in the framework of classical kinetic theory. 

For a strong high-frequency laser field a linearization with respect to the field is not possible. 

Instead, we follow the idea of Silin and assume that the influence of the collisions may be 
considered as a small perturbation compared to the external field. We decompose the current 
according to j = jo + ji and assume 

j~j° = E^ rdt'E(t'). (25) 

For a harmonic electric field, E = Eocoswt, the exponential factor in Eq. (17) can be expanded 
into a Fourier series. The current balance is given then by 

ij,, - „4 = Ke / ^ ^ E EC-^)-'.^ (3^) (^) 

OO 

xe'^' j ^ [5ee(q;a;-na;)A4(q;a;)+£f,(q;a;-na;)<S,,(q;a;)] (26) 

— oo 

with the one-component structure factors and response functions and respectively 
[8, 9]. We will assume in the following of this section that the subsystems are in local 




Vo/va,=2.00 



thermodynamic equilibrium with temperatures Tg and Tj, respectively (the influence of non- 
Maxwellian distribution functions was considered in [27], and a numerical solution of a kinetic 
equation for a strong laser field was performed in [28]). J; is the Bessel function of lih order. 
In the above equation, electron and ion functions contribute formally equal to the screening. 
The ion functions, however, are localized in the low- frequency region, i.e., for a high-frequency 
electric field, (D can be neglected in comparison with nuj. In this case, the first term in the 
brackets in Eq. (26) vanishes because J dui£^{q;LD) = 0, and for the current it follows that 

7-00 rrie ^' J i27rh)^ me mhw ''^^^^ ' 

xJn (^^) Jn-m (^^) >Cfe(q; -uu) ^^^(q) (27) 
where screening by the ions is accounted for by the static structure factor 5ji(q) defined by 

Siiici) = diDSu{ci,iD) = l + nijd^r [gu{r) - 1] e"!'!-'- , (28) 

where gu is the pair distribution function. is the exact density response function of the 
electron subsystem. Appropriate approximations can be expressed via local field corrections [9] . 

Equation (27) is clearly an expansion of the current in terms of higher harmonics of the laser 
frequency 

/t 2 °° 

dt!^B{t)= j„^(^)e'--^ (29) 

-°° m=-oo 

The Fourier coefficients jm of the current (with jm = jlm) can be easily identified from (27). 
One can show that only the odd harmonics arc allowed due to the symmetry of the interaction, 
cf. [22]. The appearance of higher harmonics in strong laser fields was first observed by Silin 
[10] and is a very interesting effect. 

Let us now consider again the cycle averaged energy (24). This quantity is now given by 

(j • E) = i r dt'm ■ E(t') = Eo • Reji = '4T^e i^eii^o) ^ . (30) 

1 Jt-T ^ ^ 

Interesting is here the restriction of the sum over m to its first term, i.e., higher harmonics do 
not contribute to the average energy exchange. 
Using the explicit expression for Reji we obtain 

/t3 oo 
^2^V^n(g)5,,(q,r,) nu;j2(^)lm£-Hq,-na;,Te). (31) 



Finally, the effective electron-ion collision frequency can be obtained using Eq. (30). 

In Eq. (31), e'^} is the full quantum mechanical dielectric function of the electrons, e.g. the 
RPA one. The physical meaning of the sum over n follows from the frequency argument nuj of 
the dielectric function. It indicates that every term of the sum describes the absorption of n 
photons (multi-photon absorption). 

In order to compare with the classical theory, it is advantageous to consider the non- 
degenerate case in which some integrations can be done analytically. For the case of a Maxwellian 
distribution function, we get 



CO 



_ 8\/27rZ^e^nerai^/mI 2 2 1 

n-i Q 



xe 



2 2 



SmekoT 



sinh 



nhjj 
2kBT 



2fesT 



(32) 



Quantum effects, marked by occur here in several places. The first place is one of the 
exponential functions describing the quantum diffraction effects at large momenta k. This 
exponential function ensures the convergence of the integral. The second place is the sinh term 
which is connected with the Bose statistics of multiple photon emission and absorption. Finally, 
quantum effects enter via the calculation of dielectric function itself. 

In the limit h ^ Q we get from Eq. (31) well-known classical results derived for the first time 
by Klimontovich [11] and later by Decker et al. [19]. The classical formulas exhibit the well- 
known problem of divergencies at large k which has to be overcome by some cut-off procedure. 
In contrast, in our quantum approach, no divergencies exist. 

Unfortunately, a further analytical simplification is only possible in limiting cases. Therefore, 
we present in the following some results of the numerical evaluation of Eq. (31) [29, 30]. In Fig. 4, 
the collision frequency is shown as a function of the quiver velocity. For comparison, results from 
Silin [10] and Decker [19] are given. Our quantum expression (31) was evaluated with the full 
dynamical and the static RPA dielectric function. The frequency dependence has been shown 
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Figure 4. Real part of the electron-ion collision 
frequency as a function of the quiver velocity 
^0 = eE/uome for a hydrogen plasma in a laser 
field {Z = l;ne = 10^^ cm-^;T = 3 ■ 10^ K; 
uj/(jjp = 5). For comparison, results of Decker 
et al. (dash-dotted line) and the asymptotic 
formulae of Silin [10] (dashed line) are given. 



already in Fig. 3. Again, we observe the resonance maximum at the plasma frequency. The 
collision frequency as a function of the coupling parameter T = (iHl)^/^ is drawn in Fig. 5. 
According to this figure our quantum mechanical results arc in good agreement with numerical 
simulations in this special situation. In contrast, the classical theories break down for F > 0.2. 
This behavior results from the cut-off procedure at large momentum k inherent in classical 
approaches. 

The influence of the Bose character of the photons [hyperbolic sine term in (32)] can be seen 
in Fig. 6. An interesting consequence is the plateau-like behavior of the collision frequency with 
a subsequent sharp drop down in dependence of the photon number. Neglecting the sinh term 
destroys this behavior. 

3. Bound free transitions. Ionization kinetics 
3.1. Basic equations 

So far we have considered transitions between scattering states. Now we take into account 
bound states and consider bound-free transitions in a partially ionized plasma, i.e., radiation 




Figure 5. Electron-ion collision frequency as 
a function of the coupling parameter T for a 
hydrogen plasma in a laser field {Z = 1; vq / Vth = 
0.2; rie = 10^^ c'ni~^;Lo/up = 3). A comparison 
is given with the theory of Decker et al. and 
with the asymptotic formula of Silin [10] (dashed 
line) is given. Furthermore, results of Cauble 
and Rozmus [31] and simulation data of Pfalzner 
and Gibbon [32] are shown. 
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Figure 6. The n— photon contributions to fgj 
vs. photon number n in a hydrogen plasma 
(ne = lO^^cm-S; Lo/up = 5; F = 0.1) for 
vo/vth = 10. Present approach (solid line), sinh 
term neglected (dashed line); classical dielectric 
theory (dotted line) 



ionization and recombination processes. Therefore, our task is to develop a kinetic description 
of these processes, especially we have to determine rate and absorption coefficients of radiation 
ionization. To investigate the problem of build-up and decay of bound states, an appropriate 
starting point of our considerations is the equation of motion of the two-particle correlation 
function 5^^, well-known as Bethe-Salpctcr equation. In order to take into account bound 
states, we have to consider the two-particle correlation function at least in binary collision 
approximation [14, 33]. Then the BSE can be written in the form [17, 3] 



7,^ 



Hah 



ab 



NabVab 



9^bM=0- 



(33) 



Here, T,^ (t) is the Hartrcc-Fock sclf-cncrgy, and Naij = 1 — — F5 the Pauli blocking factor. It 
is convenient to introduce relative and center-of-mass coordinates p and P, respectively, which 
yields 

p2 p2 _ , /etP , erp\ 



H 



ah 



2M 2/1 • 



where M and ^ are the total and reduced masses, and et = ea + the total and the reduced 
charge defined by Me^ = rribea — maeb, respectively. 



Now we take the time-diagonal limit of the two-particle correlation function, Fab{t) = 
{iKfg^^(t,t). Prom Eq. (33) there follows that [7] 



d 



ih-Fab{t) - [H°^^ + Vab,Fabit)] - [Uabit) , Fab{t)] 



[{Fait) + Fb{t)) Vab, Fabit)] - [(SHF(t) + SHF(t)) , F^bit)] = 0, 



(35) 



which describes the kinetics of the two-particle properties like bound and scattering states in the 
plasma. In Eq. (35), the first two terms represent the usual binary collision approximation, the 
third term contains the direct external field contribution, whereas the further terms arise from 
many-particle effects: Pauli blocking and Hartree-Fock (mean field -|- exchange) self-energy. 

For the further investigation, we choose, for the two-particle density matrix, the basis 
representation (9) in terms of the bound and scattering eigenstates of the field-free two-particle 
Hamiltonian defined by 



I Pat) = 



e n 



Pa) 



(36) 
(37) 



Here, a denotes a complete set of quantum numbers n, l,m in the case of bound states and the 
relative momentum p in the case of scattering states, respectively, 

_ J {n,l,m} for bound states 

p for scattering states. 



In the representation (36), we get for the equation (35) [34] 



rpaa 



(P;t) = ^ [hn'^^{P;t)FX'(P;t) - FX(P;t)hn^-'(P;t)] , (38) 



with the matrix elements given by 

(aP| Fa,{t) |pV) = FX'(P;t){27rhf5(P - P'), 

(aP| Uabit) + SH^(t) + Sr(t) ± [Fait) + Fbit)] Vab |P'a'> 

= ^J^g^"'(P;t)(27rn)^(5(P - P'), 

(aP| Uabit) + S;?F(t) + SHF(i) ± Vab [Fa{t) + Fb{t)] |PV> 

= hn'^'''(P;t){2Trh)^6(P - F'). 



(39) 



(40) 



(41) 



Here, h^l^i is an abbreviation which can be interpreted as a generalized Rabi energy, i.e., the 
field contribution renormalized by Hartree-Fock and Pauli blocking terms. Notice that Uab in 
Eqs. (40,41) does not contain the term quadratic in A (cf. Eq. (34)), because it cancels out in 
the homogeneous case in the commutator. The quantities defined in (40,41) are related to each 

other by [^Og^"'(P;t)j * = fih^'^(P;t). 

Thus, the representation chosen above leads to a transformation of (35) into a matrix 
equation. Due to the inclusion of the scattering continuum. Fab is, in principle, a matrix of 
infinite rank. If there are TV — 1 bound states, we obtain N'^ coupled equations. Even if 

one makes use of the symmetry relation FX'(P;t) = F^'"(P;t), there remain ^(-^ + 1) 
equations. 



Obviously, the diagonal elements F^^ are related to the occupation numbers of the respective 
state. Due to the external field wc have nonvanishing nondiagonal elements F^^ . They are 
connected with transitions between two states a and a' . Here, we have to consider several 
situations: 

(i) If both a and a! denote bound states, the equations describe excitation and deexcitation 
processes of atoms in the plasma. Then we recover the familiar atomic Bloch equations 
[35, 36, 37]. 

(ii) If both a and a' denote continuum (scattering) states, processes like (inverse) 
bremsstrahlung are described. 

(iii) In the case most interesting for our investigations, however, a denotes a bound state 
and a' a scattering state (or vice versa). The transitions between them, i.e. ionization and 
recombination processes, are thus included in Eq. (38). 



3.2. Plasma Bloch equations 

3.2.1. Derivation of the plasma Bloch equations Since we are especially interested in the 
ionization kinetics, in the following, we do not include (de) excitation and (de) acceleration 
processes, i.e., transitions between discrete states and transitions within the scattering 
continuum. Therefore, we consider a system with one single bound state and denote 



a 



b "bound" bound state 
/ "free" scattering state. 



We get a system of three coupled equations [34] : 



d 



ihf^F,,{V-t) = 2i j j^luv[hn^i{Y>,p-t)Ff,{V,p-t)] , 



(42) 



d 



n--{E{p)-E{p')) 



Fff{P,p,p';t) 



I 



dp 



= hnQ{P, p; i)F6/(P, p'; t) - Ffb{P, p; t)hn''^(P, p'; t), 



(43) 



^^^-(^^(P)-^^) Ff,{F,p;t)- J j^hn{l{P,p,p;t)Ff,{P,p;t) 



hn{^{P,p;t)FbbiP;t) 



dp 



:Fff{P,p,p;t)hnfJ;iP,p;t), 



Fbf{P,p;t) = [Ffb{P,p;t)]*, 



(44) 



(45) 



where E), denotes the binding energy renormalized by Hartree-Fock and Pauli blocking 
contributions. 



Eb = Eb + Jj^\Mp)\'{^f''{F,p;t) + ^r{F,p;t) 



^ [E{p) - Eb\ [Fi(P, p; t) + F2(P, p; t)\ } 



(46) 



with Fi/2(P,p;t) = Fy^{^P±p-t) and T.^J^{P,p-t) = T.^J^(^^P ±p-t). Here, 
<^6(p) = (pI'^) is the Fourier transform of the bound state wave function. 



Keeping in mind the relation (45), we have obtained a system of three equations for three 
unknown functions the physical meaning of which is quite obvious. is, in principle, equivalent 
to the distribution function of the atoms in the plasma. Fff is the binary distribution 
of the unbound electron-ion pairs. From the system (42-45) we see that the dynamics of 
both distributions is driven by the function Ff^ describing the transition between bound and 
continuum states. The latter quantity is closely connected with the polarization function. Its 
time evolution, in turn, is determined by the distributions. Notice that, even if we neglect three- 
particle collisions, ionization and recombination takes place due to the presence of the external 
field contained in the Rabi energies /iJIr. 

3.2.2. Field terms. Renormalized Rabi energies As we have seen before, the matrix elements 
of the field contribution to the Hamiltonian Uab can be represented in the form of Rabi energies 
hU^. Since we consider electron-ion pairs, Uab from Eq. (34) simplifies to 



with ipp+{p) being the Fourier transform of the scattering wave function. These relations show 
that the bare field terms are renormalized by many-particle effects in form of Hartree-Fock 
energies and Pauli blocking contributions. 

3.3. Analysis of the transition function 

In order to investigate the full ionization kinetics, the system (42-45) has to be solved 
sclfconsistently. This is, of course, a very complicated task especially due to the momentum 
dependences of the quantities which is, up to now, numerically not feasible. Therefore, we 
approach the problem by analyzing limiting cases which simplify the system (42-45) and allow 
for an insight into the underlying physics. 

Let us assume a monochromatic external field given by Eq. (2) and look at Eq. (44). 
Introducing the approximations (i) neglect of Hartree-Fock and Pauli blocking renormalizations, 
(ii) replacement of scattering states |p-|-) by free momentum states |p), and (iii) neglect of binary 
correlations in the scattering state, i.e.. 



Uab 




Thus, the Rabi energies are given by 





X (Sf F(P, p; t) + SHF(P, p; t) T [^(p) - ^(p)] [Fi (P, p; t) + F^iP, p; t)] } (48) 



^//(P>p>p'; 



t)«Fi(P,p; 



t)F2{P,p; 



t){2Trnf6ip - p'), 



this equation can be written as 



ih- - {E{p) - Eb) + ^A(t) • p F/6(P, p; t) 
= --Ait) ■ p <pbip){Fbb{'P; t) - Fi(P, p; t)F2iP, p; t)}. 



(49) 



The formal solution of Eq. (49) is given by 



ih fi 



I 

to 



X {F66(P; t) - Fi(P, p; t)F2(P, p; t)}. 



(50) 



For monochromatic fields (Eqs. (2)), the integral in the second exponent on the r.h.s. of (50) 
can be solved: 



6 I — — Vq * p 

— p • / dtA{t) = (cosa;f — cosui) . 

t 



(51) 



Keeping in mind that the field is switched on adiabatically, i.e. at to —oo, we have to 
write for the vector potential 



A(t) = lim An — 



gi(u)— ie)* g— j(aj+ie)t 



Assuming that the time dependence of the field amplitude and the distribution functions is 
weak compared to the very fast field oscillations, they can be taken out of the integral in (50). 
Then, using the relation 



^±iz cos uit 



= J2 i±^ruz)e 



we obtain 



-. oo 



ikwt 



k=—oo 



ie)t g— i(a;+je)t 



n=— oo r 

to 

xAo -pipbiP) [Fbb(P;t) - Fi(P,p;i)F2(P,p;t)] . (52) 
The time integration can be carried out, and after an index shift A; = n + Z we obtain the result 



oo oo 



2fi 

n=—oo l=—<x> 
J{{l+l)LU-ie)t 



VQ • P 

^i({l—l)ijj—ie)t 



Eb - E{p) + {n- l)huj + the Eb - -E(p) + (n + l)hjO + ihe 
X Ao • p ipb{p) [F66(P; t) - Fi(P, p; 0F2(P, p; t)] . 



(53) 



This formula can be rewritten once more by appropriate index shifts and using the Bessel 
function relations 



2n 

J-n{z) = {-lYJn{z) ; Jn-l{z) + Jn+l{z) = —Jn{z) 



(54) 



We arrive at 



^ n=-oo l=-oo 



Eb — E{p) + nhuj + ihe 



Ao • p </^6(p) [Fbb{P; t) - Fi(P, p; t)F2(P, p; t)] . (55) 



As can be seen from this result, the field causes interesting physical effects. First, the sum 
over I indicates the generation of higher harmonics of the field frequency, cf. the term e«(''^-«<:)*. 
On the other hand, the sum over n reflects the absorption or emission of multiple photons, i.e. 
multiphoton ionization, which finds its expression in the denominator E^ — E{p) + nhuj + ihe. 
Similar effects have been found, e.g., in [6, 12, 38]. 

3.4- Absorption and ionization 

3.4-1 ■ Electrical current and absorption coefficient Now we are able to consider the current 
density and other interesting physical quantities connected to j like polarization, absorption and 
ionization coefficients. For this purpose, let us come back to the balance equation (10). 
With the completeness relation for the eigenstates \Pa) and the Schrodinger equation 

1 = |P«) (P«l > {Hab - Eap) |P«) = -Vab |Pa) , 

Pa 

we express the right hand side (first line) of the balance equation (10) for the current density 
in terms of the (P; 0- follows 

diit) ,,2-p_ V /• dp f dP ^fetP erp\ 

^"'"^'^-"m^y (Wi (WtS'V^ ~/ (^pP-^ap) 

[^«(p)<^t(p)FX(P;i) - vl{p)<p^ip)FX{P;t)] =Ylab{t) . (56) 

ab 

This equation shows that the different transition processes discussed in Sec. 3.1 contribute 
to the current in a partially ionized plasma. In the following, we consider a hydrogen plasma 
as a simple model case. We should notice, however, that the theory can be applied to more 
complicated systems, too. We consider here the electron— proton part and especially the bound- 
free contribution (Q; = n = l,Q; = p and vice versa) 

Then we introduce the expression (55) into (57) and ignore the weak time dependence of 
the distribution functions and the field amplitude. Finally we integrate Eq. (56) over time and 
obtain 

oo 

m=h(.t)+ Yl j;Me^*'"* (58) 

l=—oo 

which is clearly the Fourier expansion of the current in terms of all harmonics of the field 
frequency u. The Fourier coefficients f'^e given by 

V f dp f dP e E{p) - Eb ,,e 



[F,,{P;t)-F^{P,p;t)F,iP,p;t)] J] {^) Jn-l {^) 



n=— oo 



(59) 



£^6 — E(j)) + n/ia; + i^e Ej, — E(j)) — nhuj — the 
t 

Furthermore, jo(t) = lo^i J dt'E(f') is the current of the colhsionless plasma. 

— oo 

As already mentioned above, the current determines several physical quantities, e.g. the 
polarization. Here, we will consider the energy transfer j • E between the field and the plasma. 
For the field E we assume the time dependence (2). We consider the dissipation of the energy 
averaged over one cycle of oscillation. Then we obtain after a simple calculation 

t 

(j-E) = i J dt'j(t)-E(t) = Eo-Reji(u;). (60) 

t-T 

Now we introduce the absorption coefficient for the bound-free transition by 

With the help of the Dirac identity, the substitution n — ^ —n in the second term and the 
Bessel function relations (54), after some algebra follows 



2 



V 27r ^ 3 f dp f dP E{p)-E, ^ 

X 4 {^J-^) -^b- nhw) [F,,(P; t) - Fe{P, p; t)Fp(P, p; t)] (62) 

The quantity ai,f{uj) describes the process of the absorption (emission) of n photons in the 
ionization (recombination) of atoms. In order to evaluate Eq. (62), the knowledge of the time 
evolution of the distribution functions of atoms and free particles is necessary. This evolution 
has, in principle, to be described by the equations (42) and (43). In order to obtain first 
results, we replace these quantities by thermodynamic equilibrium functions, i.e., Maxwellian 
distributions. This assumption represents a significant simplification, in particular, the influence 
of the field on the bound states is neglected. An improved determination of Ffjf, from Eq. (42) 
leads to a Stark shift of levels quadratic in the A • p interaction as was discussed also by Milonni 
and Ackerhalt [42]. 

The absorption coefficient abf{u}) is shown in Fig. 7 in dependence on the laser frequency. An 
interesting feature is the occurence of additional maxima for photon energies Hco smaller than 
the ionization energy. They are caused by two- and more-photon processes. 

If we restrict the sum in formula (62) to n = 1, our approach reproduces the relations of 
linear response theory. 



3.4-2. Rate coefficients The absorption of photons is, under the condition E{p) — E^ = nhw, 
connected with a bound-free transition, i.e. the ionization of the atoms. The process of 
ionization, of course, changes the plasma composition. From the macroscopic point of view, 
the time dependence of the densities of the plasma particles is given by the rate equation (1). If 
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Figure 7. Absorption coefficient a^f vs. laser frequency for two different temperatures. 



we restrict ourselves to radiative processes, the change of the atomic density n^(t) is given by 
(in the following, the superscript R will be dropped) 

^nA(t) = -a{t)nA{t) + l3{t)neit)np{t). (63) 

Here, ne{t) and np{t) are the densities of free particles (electrons and protons), and the 
coefficients a{t) and (3{t) arc the radiative rate coefficients. 

From the microscopic point of view, the rate equation (63) follows from the equation for the 
occupation of the atoms F{,{, (42) , 

d 2 f dv 

— Ffe6(P;t) = -sinut J ^^-|^vo •p^(,(p)ImF;b(P,p;i). (64) 
Inserting the solution (55) and using the Dirac identity, we arrive at 



^F..(P;*) = -^sm 



^ ' n=-oo/=-(x) 



X < sm 



TTCOS 



[ujt + I) <5(-E(p) -E^- nhuj) 



E}j — E{p) + nhuj 

X [Fbb(P; t) - Fe(P, p; t)Fp(P, p; t)] . ' (65) 

In the following, we assume Maxwellian distributions for the free particles, i.e. 



with the thermal wavelength Ag/p = y^^^^- Integrating Eq. (65) over P, we obtain the rate 
equation (63) for the density of the atoms nA{t). Obviously, the ionization coefficient a{t) is 



then given by the microscopic expression 

dp 



2 f dp CO CO ^ p ^ p 

a(i) = -sina;ty^vo-p|n(p)|^ ^ E Jr. {^-^) Jn+l [^-^) 

(ut + I)] d{E{p) -Eh-nhio)\, (66) 



X < sm 



l{ujt+- 



p 



Efj — E{p) + nhw 



n=—oo l=—oo 
TT COS I 



and the recombination coefficient follows from 
A? . . f dp 



n=— cxD l=—oo 



x< sm 



— E{p) + nhu 



TT cos I (^ujt + ^ j S{E{p) — Ef, — ntuo) 



xe 2^87- (67) 
with Ar = \J~^^- It is useful to average again over one period of the field oscillation. Using the 

T T 

relations ^ J dt sin cut cos sin cut sin /wt = o'^Uli) S^t the averaged coefficients 



2tt 
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"^"^ V6(P)|' E ("^)' -^n (^) S{E{P) -E,- nhw), (68) 
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Figure 8. Ionization coefficient a vs. laser frequency. 

This strongly simplified expression for the ionization coefficient is essentially equivalent to 
results from atomic physics [39, 40, 41, 42]. However, in contrast to the results of Faisal and 
Reiss, no effect of the A'^ term occurs. Effects of higher order in A following from the p • A 
coupling arc contained, in principle, in the basic equations (42-45) as mentioned already above 
in the discussion of the absorption coefficient. 

In Fig. 8 we show the ionization coefficient a as a function of the photon energy. Similar to the 
absorption coefficient we again observe a finite ionization probability for photon energies smaller 
than the binding energy due to multi-photon processes. However, we want to remark that this 
multi-photon interpretation of the ionization process is limited by the Keldysh frequency 



where / is the ionization energy. Only for frequencies larger than cok the multi-photon 
interpretation is meaningful. For frequencies smaller than lok the ionization has to be interpreted 
by the usual tunnel effect. 
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